High Energy Neutrino Flashes from Far-Ultraviolet and X-ray Flares in Gamma-Ray Bursts 
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The recent observations of bright optical and x-ray fiares by the Swift sateUite suggest these are 
produced by the late activities of the central engine. We study the neutrino emission from far- 
ultraviolet and x-ray flares under the late internal shock model. We show that the efficiency of pion 
production in the highest energy is comparable to or higher than the unity, and the contribution 
from such neutrino flashes to a diffuse very high energy neutrino background can be larger than that 
of prompt bursts if the total baryonic energy input into flares is comparable to the radiated energy 
of prompt bursts. These signals may be detected by IceCube and are very important because they 
have possibilities to probe the nature of flares (the baryon loading, the photon field, the magnetic 
field and so on). 
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High energy neutrino emission from gamma-ray bursts 
(GRBs) has been expected in the context of the standard 
internal-external scenario of GRBs. Especially, since the 
prediction of Waxman and Bahcall Jj, neutrino bursts 
in the internal shock model have been studied by several 
authors and neutrino afterglows in the external shock 
model have also been discussed 0, y| • 

The standard model of GRBs has succeeded in explain- 
ing of many observations, but there were a few outstand- 
ing questions in the study of GRBs before the launch of 
the Swift satellite (see reviews, e.g., Q)- The Swift satel- 
lite, which is an ideal mission to answer these questions, 
has presented indeed very fruitful results during the first 
several months of its operation. Especially in the x-ray 
band, Swift x-ray telescope detected an x-ray afterglow 
for essentially every burst, which showed the surprising 
behaviors that are not straightforwardly expected in the 
pre-Swift era Q. The early afterglow light curve has 
several surprising features such as an early steep decay, 
a follow-up shallower-than-normal decay, one or more x- 
ray flares and so on (see e.g., j^j). Especially, obser- 
vations show that many bursts have large x-ray flares 
superimposed on the underlying afterglow. In the early 
afterglows of XRF 050406 and GRB 050502b, x-ray tele- 
scope detected mysterious strong x-ray flares, and some 
flares such as GRB 050607 and GRB 050904 have mul- 
tiple flares ;6]. These observational results suggest the 
existence of additional emission in the early afterglow 
phase besides the conventional forward shock emission. 

In this Letter we calculate high energy neutrino emis- 
sion from far-ultravioet (FUV) and x-ray flares under the 
late internal shock model j^] ■ Our method of calculation 
using GEANT4 [g is the same as in Murase & Nagataki 
13, but quantitatively improved 93|. Now, large neutrino 
detectors such as IceCube, ANTARES, and NESTOR are 
being constructed [l3- In the near future, these detec- 
tors may detect high energy neutrino signals correlated 
with flares. 

The model. The flares typically happen hundreds of 
seconds after the trigger of prompt emissions or earlier. 



In some cases, they occur around a day after the burst. 
The observed typical time scale is larger than that of 
the prompt emission, which is St ^ (10 — 10^) s |^ |^{. 
The amplitudes of the flares are usually larger than the 
underlying afterglow component by a factor of several, 
but can be much larger. These can be even comparable 
to or higher than the prompt burst component which is 
typically Lg^^ - (10"*^ - 10^^) ergs/s around (1 - 10) 
keV band. Hereafter we take £max = (10^^ — 10^*^) ergs/s 
as a peak luminosity of far-ultraviolet and x-ray flares. 
Although some flares allow for the possibility of external 
shock processes, the general features of the flares suggest 
that this phenomenon is best interpreted as a late inter- 
nal central engine activity [5J- The variability of some 
GRB afterglows implies that the engine may last much 
longer than the duration of the bursts [iJl, although a 
possible mechanism for reactivity of the central engine is 
unknown [l3| . We suppose the late internal shock occurs 
a few minutes after the prompt 7-ray emission, powering 
a new unsteady relativistic outflow. Falcone et al. show 
that the case of GRB 050502b implies the late outflow 
has the smaller Lorentz factor, F < 20 [g- Hence, we 
assume that the typical Lorentz factors of the ejected 
material are smaller than the prompt emission, setting 
Fg ~ 10 and Ff ~ 100 as the typical Lorentz factors of 
the slow and fast shells, respectively. We can estimate 
the Lorentz factor of the merged shell F « \/FfFs ~ 30, 
and the Lorentz factor of the internal shocks can be 
estimated by F^h « {^jTi/T^ + ^VjTf)/2 - a few. This 
regenerated internal shocks are called as the late internal 
shocks W\. The typical collision radius is expressed by 
commonly used relation, r « IQ^^-''>{T /2,Qf[6t/&{l + z) s] 
cm. Of course, it should be smaller than the deceleration 
radius. Actually the exact radiation mechanism produc- 
ing the flares is unclear. Here, however, let us assume 
that flares are produced by the synchrotron emission to 
see typical parameters,. The minimal Lorentz factor of 
electrons is estimated by ^e,m ~ ee('TT'p/"^e)(rsh ~ !)• 
Since we can estimate the intensity of magnetic field by 



S = 7.7 X 10^Ge]/^_,{T,^{r,^-l)/2f^Lli;^,,rj,'r^l„ 
the observed break energy is, E'^ ~ h^yl^TeB/rrieC ~ 



0.1 keVe^ e^ 
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Here, Lm is the outflow luminosity. Therefore, the 
typical emitted energy is in the soft x-ray band. 

For numerical calculations, we simply set the break 
energy and adopt the power-law spectrum similarly to 
that of the case of the prompt emission. Although most 
of flares are actually well fitted by a Band function or 
cutoff power-law model [g, such a treatment does not 
change our results so much. We use the following ex- 
pression in the comoving frame, dn/de = nb{e/e^) for 

£""'" < e < e^ or dn/de = nb{e/e''y'^ for s'' < s < e'"'''^, 
where we set e™™ = 0.1 eV because the synchrotron self- 
absorption will be crucial below this energy and e'^^^ = l 
MeV because the pair absorption will be crucial above 
this energy |l3|. Actually we do not know the peak en- 
ergy and the lower spectral index of many flares , 7| . So 
we allow for the existence of far-ultraviolet-ray (FUV- 
ray) flares and take s^ = (0.01 — 0.1) keV in the comoving 
frame. We assume a = 1 and set /? = 2.2 j^] similarly 
to the prompt emission. The photon energy density is 
U^ = J de edn/de. The magnetic energy density and 
the nonthermal proton energy density are expressed by 
Ub = £,bU^ and Up = ^accU-y, respectively. The non- 
thermal baryon loading factor ^acc can be expressed by 
face ~ CplO(0-l/£e)j where Cp is the proton acceleration 
efficiency. Roughly speaking, face = 10 and f b = 1 cor- 
respond to €B = 0.1 if Cp ~ 1- We take r - lo"-5-i6 ^m 
with r ~ (10 — 50). Although we have done wide param- 
eter surveys, we will hereafter set the width of the shell 
to r/2r2 = 4.5 X 10"~i2 cm according to St = (30-300) 
s at z = 1, and show the two cases of r = lO^**'^ cm with 
r = 30 and r = 10^^'^ cm with T = 15. 

To obtain the pion production rate and estimate the 
maximal energy of accelerated protons, we need to take 
into account following various cooling time scales. We 
consider the synchrotron cooling time written by tgyn = 
3TOpC^/4crTTOe^pf/_B, the inverse-Compton (IC) cooHng 
time which is given by Jones [lj|, the adiabatic cooling 
time comparable to the dynamical time, and the pho- 
tomeson cooling time which is evaluated by. 
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where e is the photon energy in the rest frame of pro- 
ton, 7p is the proton's Lorentz factor, Kp is the inelas- 
ticity of proton, and eth ~ 145 MeV is the threshold 
photon energy for photomeson production. From vari- 
ous time scales, we can estimate the total cooling time 



scale by t^ _ „p^ 



tn'v + ''svn "T hn + ta 



We believe that 



not only electrons but also protons can be accelerated by 
the first-order Fermi acceleration mechanism and assume 
dup/dep oc £:~^. By the condition iace < tp, we can esti- 
mate the maximal energy of accelerated protons. Here, 
we take the acceleration time scale by iacc ~ ep/eBc. We 
also set the minimal energy of protons to 10 GeV because 



this would be around ~ TshiTT-pC^- One of numerical re- 
sults is shown in Fig. 1. 

Let us evaluate /p-y = tdyn/tp-y by A-resonance approxi- 
mation. After performing the second integral in Eq. ^, 
we can approximate by. 
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, (2) 



where cta ~ 5 x 10~^^cm^, K,p{e/^) ^ 0.2, and Sa ~ 0.3 
GeV, and Ae ~ 0.2 GeV [l|. Here, we have included 
the effect of multi-pion production and high inelasticity 
which is moderately important, and multiplied by a fac- 
tor of ~ (2 — 3) in Eq. ^ ^. Hence, we can obtain 
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where Et, ~ .be A'm.pC^T'^ / E^ is the proton break energy. 
From Eq. Q, we can conclude that a significant fraction 
of high energy accelerated protons cannot escape from 
the source without photomeson productions. In the case 
of bright x-ray flares whose luminosity is larger than 10 
ergs/s and FUV-ray flares, almost all protons accelerated 
to the very high energy region will be depleted. 

Neutrino spectrum and flux. As in the case of the 
prompt emission, we can expect high energy neutrino 
flashes from one GRB event only if the flare is nearby 
or energetic. In Fig. 2, we show an example of the ob- 
served neutrino flux from the source at z — 0.1. The 
expected muon events for above TeV energy neutrinos 
are iV^ — 0.02 events in the case of x-ray Flare (A) in 
Fig. 2. We can expect N^ = 1.4 events when the flare 
is energetic and more nonthermally baryonic [FUV-ray 
Flare (B)]. Note that, in the case of an energetic flare, 
the maximal neutrino energy will become small because 
the highest energy protons suffer from the photomeson 
cooling very much. 

We can estimate a diffuse neutrino background from 
FUV/x-ray flares for specific parameter sets under the 
standard ACDM cosmology (f^,n = 0.3, t^A = 0.7; Ho = 




log(Ep [GeV]) 

FIG. 1: Various cooling time scales and the acceleration time 
scale for FUV-ray flares with Lmax = 10^^ ergs/s, gg = 1, and 
r = lO^'*'^ cm with V = 30. Note that energy scale is measured in 
the shell comoving frame. 




FIG. 2: The observed muon-neutrino {Uf^ + U/^) spectra for one 
GRB event at 2 = 0.1. X-ray Flare (A): Lmax = lO"" ergs/s, 
^g = 1, ^^j.^ = 10, and r = IQi^-S cm with F = 30. X-ray Flare 
(B): Lmax = 10*8ergs/s, ^b = 1, €acc = 10, and r = IQi^S cm 
with F = 15. FUV-ray Flare (A): L^ax = 10*** ergs/s, ^b = 1, 
Cacc = 10, and r = IQlS-S cm with F = 15. FUV-ray Flare (B): 
Ljuax = IQSOergs/s, §3 = 0.1, Cacc = 30, and r = 10"-9 cm with 
F = 30. 



71 kms^^ Mpc^-'^). Assuming that the long GRB rate 
traces the starformation rate (SFR), we shall use the 
SF3 model of Porciani and Madau combined with the 
normalization of overall GRB rates obtained by Guetta 
et al. [13 ■ The choice of Zmax which should be larger 
than 6.3 does not affect our results so much, so we take 
^max — 11, which is the epoch when the reionization is 
likely to occur. X-ray flares are detected in at least (1/3 
- 1/2) of Swift GRBs ;5j. So we can expect that flares 
are common in GRBs. Here, introducing the ratio of 
the energy emitted by flares to that of the prompt emis- 
sion, /f = -Eflare/-£'GRB, let US estimate the neutrino flux 
analytically. First, we can express the total number spec- 
trum of accelerated protons, using ep.max ~ 10^ GeV, 
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are shown in Fig. 3 for the cases of /p^acc = (0.5 — 5). Ex- 
pected muon events from above TeV neutrinos are also 
shown in Fig. 3. Although the expected muon events 
of above TeV neutrinos from flares will be smaller than 
those from prompt bursts (shown in Fig. 3), they can be 
comparable with or exceed those from prompt bursts in 
the very high energy region above a few PeV range. 

Implications and discussions. In the near future, high 
energy neutrino signals from flares may be detected by 
IceCube and/or Auger which can provide information on 
the nature of flares. We expect such neutrino flashes 
from flares should be in coincidence with the early after- 
glow phase. Especially, some of events will be correlated 
with observed flares. These signals may be expected from 
not only long GRBs but also short GRBs, that also may 
accompany flares |^, if such flares are baryonic origins. 
However, several authors recently discussed flares may 
not be of baryonic origins but of magnetic origins '16|. 
If the outflow is much magnetized, formed shocks would 
be greatly weaken and neutrino emission would be sup- 
pressed [i^. The detection of high energy neutrinos is 
one of the tests for the origin of flares. In addition, these 
signals would give us information not only on the mag- 
netic field but also on the photon field. As demonstrated 
in FUV-ray Flare (B) in Fig. 2, too copious photon field 
will reduce the maximum proton energy and the following 
maximum neutrino energy. Moreover, such neutrino de- 
tection may include signals from FUV-ray flares, that are 
not seen due to absorption by neutral hydrogen both in 
host galaxy and in our Galaxy. Combined with gamma- 
ray large area space telescope (GLAST) mission which 
may detect sub-GeV flashes by IC Q or tt" decay, these 
high energy neutrinos may be important as a probe of 
FUV-ray emissions. 

We have estimated the flux of neutrinos from flares by 
normalizing the proton flux with the typical prompt ra- 



where fb is the beaming factor and N is the number 
of flares. Since we have already estimated fp-y in Eq. 
©, by replacing 2.hN ff,L,^^^r /T^c with the total emit- 
ted energy from flares -Efiarc, we can estimate the diffuse 
neutrino background flux from flares as follows, 
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where fz is the correction factor for the possible contribu- 
tion from high redshift sources and i?GRB(0) is the overall 
GRB rate at = where the geometrically correction is 
taken into account. /rCacc expresses the ratio of the non- 
thermal baryon energy of flares to the prompt radiated 
energy. If the nonthermal baryon energy is comparable 
to Eqy(b , we can expect a significant contribution to the 
neutrino background from flares. Our numerical results 
obtained by the same method as in our previous paper 
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FIG. 3: The neutrino background from flares. X-ray Flare; (the 
upper dashed line): Lmax = 10''® ergs/s, is = 1, /p^acc = 1, and 
r = 10^^'^ cm with F = 15; Np, = 3.7 events/yr. (The lower dashed 
line): L^ax = 10*'' ergs/s, ^b = 1, /p^acc = 0.5, and r = lO^*-® 
cm with F = 30; A^^i = 0.8 events/yr. FUV-ray Flare; (the upper 
dotted line): Lmax = 10*® ergs/s, ^s = 0.1, /pjacc = 5, and r = 



10"-® cm with F = 30; 7V„ 



1.8 events/yr. (The lower dotted 



Une): Lmax = 10*«ergs/s, ^b = 1, /p^acc = 1, and r = lO^^^-^ 
cm with F = 15; N/j, = 1.4 events/yr. GRB: taken from Q] with 
H'sh = 2xl05iergs, ^b = 1, ?acc = 10, and r = (10"_loi*-5) cm; 
Nfi = 21 events/yr. WB: Waxman-Bahcall bounds [y. 



diated energy. For the capability of detections, outflows 
may need to be largely baryon loaded. Many flares are 
likely to be interpreted as the late activity of the cen- 
tral engine. Such an energy injection is also one of the 
common interpretations for early flattening of the x-ray 
afterglow and its energy will be comparable to that of 
the prompt burst 5]. These late activities might sup- 
ply the extra nonthermal baryons. However, too large 
nonthermal baryon loading will not be plausible [2] and 
would also be constrained by high energy 7 rays. The 
high energy 7 rays could cascade in the source and/or 
in infrared and microwave background where the delayed 
emission would occur |l7l | . The delayed emission, if it oc- 
curs, would be expected in the GeV-TeV region and ex- 
tended to the keV-MeV region. In addition, there would 
also be contributions of synchrotron radiation compo- 
nents from charged particles such as pions. Although the 
detailed calculation is needed to obtain reliable spectra, 
such emission could not be detected unless the source is 
very nearby and/or the flare is energetic. [For example, 
in the case of the set with Lmax = 10''^ ergs/s, r = lO^'''^ 
cm, and F = 30, the inferred spectra [i3| imply that the 
detection by burst alert telescope would be difficult at 
z <: 0.1.] GLAST would be able to test the existence of 
high energy emission. 

At present, we do not know the total radiated en- 
ergy from flares very well. The fluences of x-ray flares 
are usually smaller than that of the bursts, but can be 
even comparable to or higher than that of the bursts. 
There are also some GRBs with multiple flares, and there 
might be a significant fraction of the FUV-ray flares. 
We also do not know the opening angle of flaring out- 
flows themselves. If the late slow outflow might have 



the larger opening angle than the prompt fast outflow, 
of which beaming factor /;, = Eq^-q / E^^°^^^ is typically 
~ (0.01 — 0.1), we can expect high energy neutrinos un- 
correlatcd with GRBs. Many parameters are still uncer- 
tain but these distributions are important for more real- 
istic predictions. For larger collision radii with r ^ 10^^ 
cm, we expect fewer neutrinos. Therefore, our evaluation 
could be maybe applied to about only a half of flares, 
and the other half would make GeV flares and very high 
energy cosmic rays. Unknown features of flares will be 
unveiled through more multiwave-length observations by 
Swift, GLAST and so on. We expect future neutrino ob- 
servations will also give us some clues on the physical 
parameters of flares. 

Waxman and Bahcall [ij predicted neutrino burst un- 
der the assumption that GRBs are the main sources of 
ultra- high-energy cosmic rays (UHECRs). In the case of 
the prompt bursts, the optical thickness for the photome- 
son production can be smaller than the unity especially 
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cm and the 



at larger radii r > 10'''(£:^"°h/10''^ ergsj 
UHECRs can be produced in such regions. In flares, it is 
more difficult to generate UHECRs especially in the case 
of FUV-ray flares where the optical thickness for the pho- 
tomeson production can be larger. So far, we have not 
taken account of neutrino oscillations. Actually, neutri- 
nos will be almost equally distributed among flavors as a 
result of vacuum neutrino oscillations [ij . 

We thank the referees, K. loka, K. Toma, K. Asano, 
and Z.G. Dai for many profitable suggestions. S.N. is 
partially supported by Grants-in-Aid for Scientific Re- 
search from the Ministry of E.C.S.S.T. of Japan through 
No. 16740134. 



[1] 
[2] 
[3] 



[4] 
[5] 

[6] 

[7] 
[8] 



E. Waxman and J. Bahcall, Phys. Rev. Lett., 78, 2292 
(1997); Phys. Rev. D, 59, 023002 (1998). 
K. Murase and S. Nagataki, Phys. Rev. D, 73, 063002 
(2006). 

E. Waxman and J. Bahcall, ApJ, 541, 707 (2000); Z.G. 
Dai and T. Lu, ApJ, 551, 249 (2001); CD. Dermer, ApJ, 
574, 65 (2002); CD. Dermer and A. Atoyan, Phys. Rev. 
Lett., 91, 071102 (2003); D. Guetta, D.W. Hooper, J. 
Alvarez-Muniz, F. Halzen, and E. Reuveni, Astropart. 
Phys., 20, 429 (2004); K. Asano, ApJ, 623. 967 (2005). 
B. Zhang and P. Meszaros, IJMP. A, 19, 2385 (2004); T. 
Piran, Rev. Mod. Phys., 76, 1143 (2005). 

astro-ph/0506157 B. Zhang 

( 2006)1 



al., 
354 



P.T. O'Brien et al. 



S. Kobayashi et 

et al, ApJ, 642 

astro-ph/0601125 

D.N. Burrows et al, Sci., 309, 1833 (2005); S. Barthelmy 

et al, Nat., 438, 994 (2005); G. Cusumano et al, Nat., 

440, 164 (2006); A.D. Falcone et al, ApJ, 641, 1010 

(2006); P. Romano et al, A&A, 450, 59 (2006). 

Y.Z. Fan and D.M. Wei, MNRAS, 364, L42 (2006); Y.Z. 

Fan and T. Piran, MNRAS, 370, L24 (2006). 

S. Agostinelli et al, Nucl. Instrum. Methods Phys. Res., 

Sect. A, 506, 250 (2003). 



[9] S. Schadmand, Eur. Phys. J. A, 18, 405 (2003); Particle 

Data Group, http://pdg.lbl. gov ^ 

[10] ANTARES Collaboration, astro- ph/9907432| (1999); 

P.K.F. Grieder et al, Nuovo Cimento, 24C, 771 (2001); 

J. Ahrens et al, Astropart. Phys., 20, 507 (2004). 
[11] K. loka, S. Kobayashi, and B. Zhang, ApJ, 631, 429 

(2005). 
[12] D. Proga and M. Begelman, ApJ, 592, 767 (2003); A. 

King et al, ApJ, 630, L113 (2005); R. Perna, P.J. Ar- 

mitage, and B. Zhang, ApJ, 636, L29 (2006). 
[13] K. Asano and F. Takahara, PASJ, 55, 433 (2003); Z. Li 

and L.M. Song, ApJ, 608, L17 (2004). 
[14] F.C Jones, Phys. Rev., 137, B1306 (1965). 
[15] C Porciani and P. Madau, ApJ, 548, 522 (2001); D. 

Guetta, T. Piran, and E. Waxman, ApJ, 619, 412 (2005). 
[16] B. Zhang and S. Kobayashi, ApJ, 628, 315 (2005), Y.Z. 

Fan et al, ApJ, 635, L129 (2005); Z.G. Dai et al., Sci., 

311, 1127 (2006); D. Proga and B. Zhang, MNRAS, 370, 

L61 (2006). 
[17] CD. Dermer and A. Atoyan, A&A, 418, L5 (2004); S. 

Razzaque, P. Meszaros, and B. Zhang, ApJ, 613, 1072 

(2004). 



